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(N ; Abstract 

q . A linear stability analysis of twisted flux-tubes (strings) in an SU(2) semilocal theory - 

an Abelian-Higgs model with two charged scalar fields with a global SU(2) symmetry - is 
carried out. Here the twist refers to a relative phase between the two complex scalars (with 

psj ■ linear dependence on, say, the z coordinate), and importantly it leads to a global current 

flowing along the the string. Such twisted strings bifurcate with the Abrikosov-Nielsen- 
Olesen (ANO) solution embedded in the semilocal theory. Our numerical investigations 
of the small fluctuation spectrum confirm previous results that twisted strings exhibit in- 
stabilities whose amplitudes grow exponentially in time. More precisely twisted strings 
with a single magnetic flux quantum admit a continuous family of unstable eigenmodes 
with harmonic z dependence, indexed by a wavenumber k G [— k m ,k m \. Carrying out a 

fS| . perturbative semi-analytic analysis of the bifurcation, it is found that the purely numer- 

ical results are very well reproduced. This way one obtains not only a good qualitative 
description of the twisted solutions themselves as well as of their instabilities, but also a 
quantitative description of the numerical results. Our semi-analytic results indicate that 
in close analogy to the known instability of the embedded ANO vortex a twisted string is 

g ; also H My to expand in size caused by the spreadtng out of its magnetic fto, 

G\ . 

Q , Keywords: stability of semilocal strings 
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Introduction 



Topological defects, such as domain walls, strings and monopoles arise in many field theory 
models with spontaneous symmetry breaking. These objects are expected to be created during 
phase transitions, and are likely to play an important role in the early universe and more 
generally in cosmology [T]. Typically these defect solutions are stable because there is an infinite 
energy barrier separating them from the vacuum. An important class of defects is constituted 
by line-defects called cosmic strings. A cosmic string is typically a flux-tube, enclosing a 
certain number of magnetic flux quanta. In the plane orthogonal to its direction a cosmic 
string corresponds to a vortex solution. Vice versa, any vortex in the plane can be extended to 
a straight string in the orthogonal direction to the plane. The well-known Abrikosov-Nielsen- 
Olesen (ANO) vortex of the Abelian Higgs model serves as a prototype straight cosmic string. 
ANO vortices are characterized by an integer winding number, n, determining their magnetic 
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flux, and by the mass ratio, (3 = ml/ ml (m s resp. m v denoting the mass of the scalar resp. 
vector fields). 

Rather interesting, cosmic string-type defect solutions have been found in semilocal theories 
where there is no topological stability [2j[3j|4]. An important class of semilocal theories is pro- 
vided by Abelian Higgs models with a suitably extended scalar sector, the simplest one being 
two complex scalar fields admitting global SU(2) symmetry. Semilocal models are quite inter- 
esting since both global and local symmetries are simultaneously broken; for a comprehensive 
review we refer to [5j. Abelian Higgs models with an extended scalar sector contain ANO- 
type vortices, which correspond to a simple embedding. Quite remarkably semilocal models 
exhibit stable vortices, despite the first homotopy group of the vacuum manifold being trivial 
[21 [3]. The linear stability of ANO vortices embedded into SU(N) symmetric semilocal models 
has been examined in Ref. [3], by computing their small fluctuation spectrum. It has been 
found that embedded ANO vortices are stable only if (3 < 1. For {3 > 1 there is a family of 
energy lowering eigenmodes, parameterized by the ^-direction wave number, k G [—k m ,k m ], 
with the lowest lying eigenmode being z-independent. This homogeneous unstable mode can 
be seen to correspond to a "magnetic spreading" instability [21 [7]. The long time dynamics of 
the instability (i.e. the numerical solution of the full nonlinear field equations) has also been 
studied, in Ref. [TO]. There it is found that for (3 > 1 embedded ANO vortices undergo indeed 
a homogeneous expansion as expected from the analysis of Ref. [3]. 

The case (3 = 1 is quite special, for here, instead of there being a unique (and stable) vortex 
solution satisfying the first order Bogomol'nyi equations [6], there is a continuous family of 
them having the same energy [3]. The Standard Model counterparts of the embedded ANO 
solutions are the Z-strings, which have been shown to be stable for sin 2 W > 0.9 [8]. 

Recently, a new family of current carrying, "twisted" string solutions have been discovered in 
the SU(2) symmetric semilocal model [9|. These solutions are characterized by a relative phase 
difference, exp(iuz), between the two components of the scalar field, where z is the coordinate 
along the string and u is the twist. Twisted strings exist only for (3 > 1 where the embedded 
ANO solutions are unstable. Quite remarkably the energy of twisted strings is lower than 
that of the embedded ones (as a matter of fact the energy per unit length of twisted strings 
is a monotonously decreasing function of the twist u>). The instability of the embedded ANO 
vortices for f3 > 1 can be seen to correspond to a bifurcation with twisted strings. Clearly the 
problem of stability of twisted strings is an important one. 

A linear stability analysis of twisted strings has been presented by Ref. [16j, and a family of 
unstable modes has been found. More precisely for all values of the twist u, the lowest eigenvalue 
belongs to a nonzero value of the wave number k of the eigenmode in the longitudinal direction. 
In particular there are no homogeneous (z-independent) negative energy eigenmodes. Noting 
the analogy between the harmonic z-dependence of the instability mode and the Plateau 
Rayleigh instability in hydrodynamics (i.e. the fragmentation of a fluid stream into droplets), it 
has been argued in Ref. [H] that the unstable eigenmode of twisted strings signals its breakup 
into small droplets. 

In this paper, we present a detailed analysis of the stability of twisted vortices. Working in a 
different gauge than Ref. [16] our work provides an independent check of the results of Ref. [16J. 
Both Ref. [16] and the present paper apply linearization analysis. In the present paper, following 
Ref. [H] . we work in the background field gauge. By solving the small fluctuation equations 
numerically for a large range of the parameters of the problem we confirm the numerical results 
of Ref. [16J, although we also find some small discrepancies which do not effect, however, the 
validity of the numerical results. We also find that twisted strings possess a family of unstable 
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modes, parametrized by their wave number in the z direction, k G [— k m , k m \. These results are 
in agreement with those of Ref. [IE]. In addition to a full fledged numerical approach, we have 
also developed a semi-analytical description of the instability eigenmodes of twisted vortices for 
values of the twist near the bifurcation point. In our view this sheds some light on the nature 
of the instability of twisted strings. The description of the bifurcation also makes it possible to 
obtain the eigenvalues and eigenmodes of embedded strings as a deformation of those of twisted 
strings (with the deformation parameter being a function of the difference of the twist u> and 
the bifurcation point u^). We note that twisted strings with periodic boundary conditions are 
also unstable in the semilocal model. 

Clearly a linear stability analysis is not sufficient to draw definite conclusions on the issue 
of the final state of twisted strings. There are nevertheless two important properties of the in- 
stabilities which can form the basis of some speculations on the long-time dynamics of twisted 
vortex instability. Firstly, not far from the bifurcation with the ANO vortices, the unstable 
modes of twisted vortices and those of the embedded ANO ones are very similar. The field com- 
ponent, which dominates the eigenmodes close to the one corresponding to the largest negative 
eigenvalue is nearly homogeneous. Secondly, independently of the absence of z-independent 
instability eigenmodes for twisted strings the z -independent instability mode of the embedded 
ANO vortices is still an energy lowering perturbation. In such sufficiently general 

initial perturbation (like a lump on the string, see Subsection 12.41) overlaps with the eigenmode 
corresponding to the most rapidly growing ones, those being almost ^-independent. Therefore 
we predict that twisted strings whose twist is close to the value at the bifurcation point start to 
expand nearly homogeneously very similarly to the instability of embedded ANO vortices. For 
twisted strings whose twist is is far from the bifurcation value this expansion may become local, 
in the sense that expanded lumps may form on the strings. The long-time dynamics, however, 
is expected to be rather different for twisted strings as compared to the embedded ANO ones. 
This difference is due to the global current flowing in twisted strings. We have found that to 
first order in perturbation theory the current remains localized to the string. Unless the current 
is completely carried away by radiative effects which come from higher orders in perturbation 
theory, in contrast to the expansion of the embedded ANO strings [TO], the expansion of twisted 
ones cannot go on indefinitely. It is also conceivable that the instability of static (or stationary) 
twisted strings signals that there are time dependent (oscillating) breather-like states. We also 
present some arguments, that twisted strings are unlikely to break up into small droplet-like 
configurations. 

The outline of the paper is as follows: in Section[Uwe recapitulate the main characteristics of 
the SU(2) semilocal theory and its vortex solutions. In Subsection [L3] we give a semi-analytical 
description of the bifurcation of twisted vortices with the embedded ANO vortex. In Section [2] 
we turn our attention to the stability problem of the twisted vortex solutions, and present the 
numerical results. In Subsections 12.21 and 12.31 the behavior of the linearized equations near the 
bifurcation is studied. In Subsection 12.41 the properties of the eigenmodes are analyzed, and 
the possible scenarios for the dynamics of the instability are presented. The technical details 
of the calculations of Section [2] are relegated to Appendix lAl In Appendix iBl a summary of the 
perturbations of embedded ANO vortices is given. 
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1 The SU{2) semilocal theory and its vortices 



1.1 The SU{2) semilocal theory 

The 4 dimensional (4D) action of the SU(2) semilocal theory can be transformed by suitable 
rescalings, rendering the rescaled fields and coordinates to be dimensionless, and setting the 
charge of the scalar particle to unity, to the form: 

S = Jd A x j-^^ + ^tM-^rt-lfJ . (1) 

Where F pu = d p A u - d u A p} = <9 M $ - L4 M $, $ T = (0i,0 2 )- The signature of the flat 

Minkowskian metric used here is (+,—,—,—). Because of the spontaneous breaking of the 
U(l) gauge symmetry the physical spectrum contains two massive particles, a scalar and a 
vector whose mass ratio is given by (3 = m^/m^. The fields transform under the U(l) gauge 
symmetry as 

A^A^ + d^x), $^e iA W$, (2) 

while the complex doublet, $ T = (</>i,</>2), transforms as the fundamental representation of the 
global SU(2) symmetry. The Euler-Lagrange equations following from the action (DQ) can be 
written as: 

d?F p , = z{(D M $)t$ - &D^} , (3a) 
D P D P $ = (3(1 - |$| 2 )$, (3b) 

where |$| 2 = $t<& = |0i| 2 + |02| 2 - Any solution of the Abelian-Higgs model with a single scalar 
field, (Ap , 0i), can be embedded into the semilocal model simply by putting (A p , <£> := 0i$o) 
where <£>o is a constant SU(2) doublet of unit norm, $J$o = 1- 

The SU(2) symmetric semilocal theory has the following conserved Noether currents 

j; = -i ((D^ a )*T^ b - , (4) 

where a = (0, a) and T a = (l,r a ) with r a being the Pauli matrices. The currents generate 
the global SU(2) while j° the local U(l) symmetry. 

We shall also consider the special theory corresponding to the (3 — > oo limit. In this limiting 
case, the scalar fields are constrained by = 1, and the limit theory is nothing but a gauged 
CP 1 -model. The constraint can be taken into account by replacing the scalar potential in the 
Lagrangian of Eq. (J3j) by a Lagrange multiplier term: 

Coo = -\f^ + {D^yD^a - A(|0| 2 - 1). (5) 



The equations of motion in this case are very similar to the regular case, 

(6) 



D p D^ a = -A0 a 



with 



A = faD^a = MD^D^aY . (7) 



4 



1.2 Vortex solutions 



To obtain twisted vortex solutions of the theory (PQ), one needs to write up the most general 
axially symmetric Ansatz [9]. Here, "axially symmetric" is meant in the general sense, that 
rotations around the symmetry axis, or translations in its direction can be compensated by a 
suitable gauge transformation [12]. 

Let us choose the z = x 3 coordinate along the symmetry axis. Choosing a suitable reference 
frame the configuration can be rendered static, which choice shall be assumed from now on. 
The axially symmetric configuration has rotational symmetry in the (x\ , x 2 ) plane, and is then 
described by the following general Ansatz: 



0i(r,0,z) 
0a(r,0,2) 
A#(r,$,z) 



/i(r)e« 

/ 2 (r)e fm V^, 

na(r), 

ua 3 (r), 



(8) 



with A = 0, A r = 0. The parameter u, which describes the z dependence of the relative phase 
of the two scalar field components, is the twist. 

In what follows we shall consider the simplest and most important class of configurations 
with m = and with n = 1, but of course our analysis can be easily extended for the general 
case. The vortex profile functions in (jHJ) obey the equations 



1 

-(ra 

r 



a)' 



2a 3 |/| 2 -2/ 2 . 



r - = 2/ 2 (a-l) + 2/ 2 a. 



r 



2)' 



fi 
h 



a) 2 n 2 



[na 



+ a/a 2 -/3(1 - 



(9) 



+ c 2 (l-a 3 ) 2 -/3(l-|/r) 



Regularity at the origin r 

j 



is ensured by the boundary conditions a 



0, A - 0, f 2 ^ c, 

a 3 — > c, where c, d are constants, while for r — > 00 we demand a — > 1, /1 — > 1 /2, a 3 — > 0. 

The embedded Abrikosov-Nielsen-Olesen (ANO) solutions correspond to / 2 = ct 3 = in 
Eqns. (JHJ). Besides these, new classes of vortex solutions have been discovered. First, Ref. [3] 
found a one-parameter family of solutions for (3 = 1 with f 2 7^ and a 3 = 0. The parameter 
of the vortices in this family of solutions can be identified with their width (denoting the 
parameter with £, the radial decay of the magnetic field is 1 — (r/^) -2 ). These solutions satisfy 
a Bogomol'nyi type energy bound [2], therefore their energy is degenerate, and they obey 
simpler, first order field equations. 

Another one-parameter class of solutions has been discovered in [9] for (3 > 1 with both / 2 
and a 3 7^ 0. These solutions are parameterized by the value of their twist, < uj < Wb(/3). At 
uj = u h these twisted vortices bifurcate with the embedded ANO vortex (for some numerical 
values of 0Jh{P) for a range of (3 see Ref. [9]). There is a similarity between this family of vortices 
and the one of Ref. [3]. The parameter u also characterizes the width of vortices: a lower value 
of u corresponds to a more diluted vortex with a slower radial decay of the magnetic field, 
and a higher value of the scalar condensate at the origin. (However, one should note, that the 
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radial decay of these solutions is different from that of the (3 = 1 ones, as here, f 2 decays with 
exp(— ur).) 

An important property of twisted vortices is that there is a global current flowing along them. 
In the coordinate system fixed by the Ansatz (jHJ), the current is the third isospin component 
of the semilocal current (J4]) in the z direction. The total current flowing along the string can 
thus be written as 

I 3 = Jrdrd&jl. (10) 

As the twist uj decreases, the current increases and the twisted vortex becomes more and more 
spread out [9]. It appears that the limit uj — > is somewhat singular, nevertheless in this limit 
twisted vortices converge pointwise to a configuration which is precisely the large width limit 
of the (3 = 1 semilocal vortices of Ref. [3]. 

Another important property of twisted vortices is that they have lower energy than the 
embedded ANO solutions. This indicates that such twisted vortices may be preferred by phys- 
ical processes over the embedded ANO vortices. In fact, for fixed (3, the energy monotonously 
decreases as uj decreases, i.e. as we go further away from the bifurcation [9]. 

In the [3 — > oo limiting theory fx and f 2 are related by the constraint / 2 + /f = 1. A 
convenient way of parameterizing them is 

/i = sin0, / 2 = cos0, (11) 

(12) 

sin(2#) . 

In this case the embedded ANO solution has infinite energy, whereas the twisted vortices are 
of finite energy and they exist for all values of uj [9]. 



which yields the profile equations 



r 



(ran 



a 



2 [o 3 - cos 2 6] 



r [ — ) = 2(a — sin 2 ■ 

r 



- W = £ 

r 2 



(2a 3 - 1)uj' 



n 



(2a - 1) 



1.3 Bifurcation of the embedded ANO solution 

It is by now well known (3] that the embedded ANO vortices are unstable to small perturbations 
of the / 2 variable. This instability corresponds to the bifurcation of the ANO vortices with the 
twisted vortex solutions [9]. The systematic expansion of a twisted vortex near the bifurcation 
point can be then written as: 

A = A (0) + e 2 f[ 2) + ... 
/ 2 = e / 2 (1) + e 2 /f + ... 

a = a (0) +e 2 a (2) + _ ( 13 ) 

a 3 = e 2 <4 2) + • • • 
uj = uj h + euJi + e 2 uj 2 + ■ ■ ■ 

where , f[ ^ denotes the ANO vortex, whose equations can be read off from equations Q 
by putting f 2 = a 3 = 0. In the above expansion, we have omitted some terms linear in e. When 
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solving the equations of the field components, these turn out to be zero (i.e. their equation 
is a homogeneous linear equation). As we shall see the only first order term is / 2 . This is 
physically plausible too: the bifurcation is parameterized by the growing of the condensate / 2 . 

The equations for the field components can be obtained easily by substituting the above 
expansion into the vortex profile equations ((9|). To first order in the e expansion the only 
non-trivial equation determines the second component of the scalar field, 

(nam 2 



(The operator should not be confused with the /i = 2 component of covariant derivative 
D^) In second order the equation for f^p can be written as 



7 ("= 



(1)' 



P (i - (A (0) ) 2 ) + 



10-, 



f: 



(i) 



0. 



(14) 



(I?J )) +^)/ a w = -2 WbWl / 2 w . (15) 
Eq. (fT5l) contains only a resonance term on the right hand side, therefore the solution of this 



r(2) 



equation is /. 



(2) 



and the removal of the resonance term yields uj\ 

4 2) T/r-2 (/ r) 2 4 2) = -2(/ 2 {1) ) 2 , 



0. 



rfi J /r 



ran 



n 2 (l-a^) 2 



,(2) /S 



- /3{1 - 3(J®)* 
2(/(°) )2 a (2) _ 4 ( a (o) 



i)A (0) A (2) 



2n 2 (l -a(°) N 



,(2) 



/?A (0) (/ 2 (1) ) 2 



2a(°)(/ 2 (1) ) 2 



The absence of the resonance term in Eq. (TT6l) gives 



uO-2 



2o; b ||/ 2 (1) | 



(1) 

2 ) 



2na(°) 



,2^(2) Al) 



(0) ,(2) ,(1) 



f(l)\3 



where the scalar product and the norm is defined as 



(f,9)--- 



rdrfg 



and 



if, f) ■ 



(16) 



(17) 



(18) 



Identifying /2(C)) with the perturbation parameter, e, at the bifurcation point, we obtain that 



LO2 



-(uj -Ub) + ... . 



(19) 



We have compared the the results of the above perturbative bifurcation analysis, for /2(0) as a 
function of uj with numerical results on Figured! A fitted curve, of the form u 2 l ^ 2 {oJb — uj) 1 / 2 + 
ei(u)b — to) , with oj 2 and €\ being the parameters fitted, as suggested by perturbation theory, 
is also shown. See also Table El Another nice comparison of the twisted and ANO vortices 
can be seen on Figures EH31 We have drawn no figures of fi and a. Their perturbation series 
starts with a zeroth order ANO term, and the next correction is of order e 2 , therefore, the 
difference between the ANO background, the second order perturbative result and the exact 
twisted vortex profile function is quite small. 

As we can see from both the comparison of the e - uj dependence and the vortex profile 
functions, the above perturbative method gives a good approximation of twisted vortices, even 
relatively far from the bifurcation, i.e. for e 



0.2. . .0.3. In Subsections 12.21 and 12.31 we will use 
this perturbatively obtained vortex background to explore the relation between the instabilities 
of the embedded ANO vortices and those of the twisted ones. 
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2 4 6 8 10 12 14 16 18 



Figure 2: The twisted vortex profile function f 2 for (3 = 1.25, uj = 0.13. For comparison, we 
have also plotted fi and a in this figure (note the scaling by 1/4 of fi and a). 

2 Linear stability analysis 
2.1 Perturbation analysis 

Here we present a stability analysis of the twisted vortex solution based on the linearization of 
the field equations about the vortex. Let 4> a = 4> a ,bg + ^4>a and = A^ ^g + e<5A M , with 4> a ,bg 
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Figure 3: The twisted vortex profile function 03 for (3 = 1.25, to = 0.13 

and A^bg being the background solution. The index "bg" will be dropped in what follows; the 
total perturbed fields will not appear later on. Here, we will summarize the main properties 
and the method of solution of the linear equations governing Sep a and SA^. Following the 
footsteps of Ref. [H] we will use the background field gauge. The main differences between the 
present case and that of the ANO vortices is that here there are more field components and the 
z-dependence does not decouple, which makes the problem computationally more demanding. 
The main advantage of the background field gauge is the simplicity of the resulting equations. 
Its main drawback is, that, unlike the method used by Ref. [IE], it does not fix the gauge 
completely, which makes an analysis of ghost modes necessary. 

The perturbation equations can be derived from a linearization of the field equations. The 
resulting equations, taking into account that the background solution satisfies the Lorentz gauge 
condition d^A tl = 0, can be cast into the form 




(20) 



where the operator D is 



'A + 2i5 ab A^d, 



-id v (f>* a + 2A v (f)* a + z< 



-/3<fib4>a 

D 1 - 2i8 ab A»d tl - 
id v <p a + 2A v <f> a - 



2id u 6 b + 2A L 



-2id^t + 2A^* h - i<t>* h d, 



<9/A + g^n 



(21) 



with Di = -5 ab D + A^5 ab - (3(\<p\ 2 - l)5 ab and □ 

In the P — > 00 limit, similarly to the regular case, the Lagrangian is expanded to second 
order in 6<j) a , S(p* and SA^. Denoting the zeroth and first order terms of D^D^(p a with \a and 
ip a , equation (16T1) can be written as ip a = 0. Now, the linearization of the Lagrange multiplier 
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has to be added. The resulting equations of motion are 

^ - (A {o) 50 a + A«0a) = 



(22) 

r a - (A (0) 5C + A«C) = o, 

where 

A (0) = \ {<t>* a Xa + <PaX* a ) , 

l (23) 

A (1) = ^ ( S €Xa + Ma + 5<t>aX*a + M*a) ■ 

The perturbation equations are invariant under infinitesimal gauge transformations of the 
form 

5(f) a -> 5(j) a + ix4> a , 
SAp^SAp + dpX- 

This gauge freedom can be dealt with by adding a gauge-fixing term ||F(A)| 2 to the Lagrangian, 

F(A) := dpSA 1 * + i{5(/)*J a - <f>* a 8<t> a ) = 0. (25) 

This helps one to get rid of the terms containing first order derivatives in D, yielding the 
perturbation matrix 

'£>i + 2i5 ab A^ -{13+ l)M* a ~{P - l)060a 2id^ b + 2A^ b \ 

-{(3 - l)0 fe *0: D x - 2i6 ab A»d lt -{(3 + iWJa -2id^* b + 2A^ h , (26) 

-2id v <f>l + 2A v <f>l 2id v (j )a + 2A v (j )a \<P\ 2 2g^ + g^U J 

where D x = -S ab D + A^5 ab - /?(|0| 2 - l)5 ab . 

It can be easily verified that the equation of SA , 

d ll & 1 5Ao + 2\<j)\ 2 5Ao = 0, (27) 

decouples, just as in the case of the ANO vortex [llj. 

The gauge condition (1251) does not fix the gauge completely, but still allows gauge transfor- 
mations whose generating function satisfies the ghost mode equation 

<9^x + 2|0| 2 x = O. (28) 

The ghost modes cancel the SA modes, whose Eq. (1271) is identical to (1281) . and a part of 
the spectrum of the gauge fixed operator fl26l) . In the case of the embedded ANO vortices, a 
straightforward analysis of ghost modes is possible, see Appendix iBl 

Let us now apply a Fourier transform to the perturbation equations (1261) in the t and z 
variables. Introducing \I/ = {5<fri,5<frl,5(t)2,8(t>2,5Ai,8A 3 ,5AQ) T . 

= / dkdQe i{nt - kz) y. (29) 



This way, the perturbation equations can be brought into the form 

= tt 2 y. (30) 
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See Appendix [A] for the details of the calculation and the matrix M.. The angle variable $ can 
be separated in a similar fashion, by the Fourier series expansion 

* = J2 *te m - (31) 

i 

This yields the desired eigenvalue problem 

M^ = n 2 ^, (32) 

where is an ordinary differential operator in the radial variable r (see equation (16T1) in 
Appendix [AT) . A similar expansion can be done for the ghost modes. 

The numerical method used for solving equations (1321) as well as the background equations 
(jHJ) was shooting to a fitting point. The resulting eigenvalues as a function of k can be seen on 
Figure 0] for (3 = 1.25, 2, 2.5 and the limit (3 — > oo. There are two unstable modes, one which is 
the deformation of the unstable mode in s 2 of the embedded ANO vortex, which has a negative 
eigenvalue for < k < k m , and another one, which is the deformation of the unstable mode in 
s 2 of the embedded ANO vortex and has a negative eigenvalue for > k > —k m . These two 
modes are related by a transformation of k — > —k and interchanging the fields si,s 2 ,a + and 
the conjugate fields s*,S2,a_. The numerical results exhibit this symmetry to a rather high 
accuracy. 

In Figure SI we have plotted the dispersion relation of the unstable modes for < k < k m . 
A closed form of this is not known, however, the approximation 

n 2 = n 2 min + n 2 2 (k-k miQ ) 2 (33) 

is empirically found to be very good. Note that fc mm , the wavenumber corresponding to the 
minimal eigenvalue, should not be confused with k m standing for the maximal value of the wave 
numbers of the unstable modes. In Figure [5] one can see f2 min and k min as a function of uj for 
(3 = 1.25. The same data, for (3 = 1.25, 2 and 2.5 is summarized in Tables CD, [2] and [3l In the 
case of the embedded ANO vortex, the form (1331) of the dispersion relation is exact. In Ref. 
[16] it is claimed that the minimum value is achieved at k = ±a>, but for smaller values of u; we 
have found a considerable deviation from /c min = ±u (see especially the curves corresponding 
to (3 — > oo). 

Close to the bifurcation, our numerical results agree with those of Ref. [16]. However, for 
smaller values of oj, there is some discrepancy which we attribute to numerical errors. 



u 


n 2 . 

mm 


n 2 


k ■ 
"'min 


0.07 


-0.00105 


0.73927 


0.05418 


0.08 


-0.00173 


0.80044 


0.06383 


0.09 


-0.00262 


0.84997 


0.07352 


0.11 


-0.00591 


0.94158 


0.09542 


0.12 


-0.00906 


0.97588 


0.10882 


0.13 


-0.01390 


0.99548 


0.12435 


0.136 


-0.01801 


0.99984 


0.13513 



Table 1: Twisted vortex perturbations: parameters of the dispersion relation for (3 = 1.25 
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0.27098 


0.32 
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0.30837 
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Table 2: Twisted vortex perturbations: parameters of the dispersion relation for (3 = 2 
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0.42 
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Table 3: Twisted vortex perturbations: parameters of the dispersion relation for (3 = 2.5 
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0.266817 
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2 


-0.402749 
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3 
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Table 4: Twisted vortex perturbations: parameters of the dispersion relation for (3 — > oo 
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1.25 


0.13694 
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2.0 
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-0.17824 


0.32989 


-0.17815 


-0.1088 


0.488 
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0.3914 
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2.5 


0.42744 


-0.21244 


0.42744 


-0.21225 


-0.1827 


0.720 


-0.1826 


0.6003 
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Table 5: Perturbative and numerical data of vortex instability modes. Here, the subscript "fit" 
denotes data obtained by fitting a parabola to the numerically obtained dispersion relation, 
close to the minimum. 
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Figure 4: The dispersion relation of the unstable mode 




Figure 5: Minima of Q 2 and its location for (3 = 2 as a function of u; 



It is important to note, that the modes discussed in this section are all part of the physical 
spectrum. This can be verified by solving the ghost equation (l28l) . The ghost modes we 
obtained do not have nodes and correspond to a positive eigenvalue, therefore there cannot be 
lower eigenvalue modes. This shows, that the instability modes presented in this paper are all 
physical (i.e. there are no negative eigenvalue ghost modes which could cancel them). 

For the (3 — > oo case the results of the numerical calculations can be seen (with the 
parametrization used in the (3 < oo case) in Table QJ See also the dispersion relation in 
Figure QJ In this case, there are no ghost mode solutions (note also that in the regular case, 
the ghost mode eigenvalues increase for larger values of (3). The (3 — > oo data shows that the 
instability modes persist for all values of (3. For (3 — » oo uj h — > oo, thus the (3 — > oo data can be 
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Figure 6: Minima of Q 2 and its location for (3 — > oo as a function of to 





Figure 7: f^ =tl) and its comparison with PT and the dependence of the minimal Q 2 on e for 
(3 = 1.25 



used to study vortices with values of uj <C uJb- Our data indicates that for (3 fixed, the family 
of unstable modes parameterized by k G [—k m , k m ] also persists for all values of < uj < uj\>. 



2.2 Bifurcation analysis of the perturbation operator 

For values of the twist uj ~ uj^ we have expanded the vortex backgrounds in a suitable parameter 
e. In this paragraph we give the second order expansion of the perturbation operator and its 
eigenvalues in the bifurcation parameter e. For details, see Appendix IA.31 
The leading order equation of the mode pertinent to our problem is 

D^s 2 = (n 2 )^s 2: (34) 

which yields the unstable mode of the embedded ANO vortex. Let us note here, that for k — uj, 
this instability mode is z-independent. For other values of k, its z-dependence is harmonic. 

Let us now look at the first order correction to the eigenvalue. This correction, as known 
from quantum mechanics is given as 

(^ 2 ) (1) = ^ (* 2 ,^V) = 0. (35) 

Therefore, to obtain an w-dependent eigenvalue, one must go to the second order in pertur- 
bation theory. Contributions quadratic in e come from two sources: from 0(e 2 ) terms in the 
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perturbation operator (i.e. corrections to the background functions of order e 2 ), and from second 
order perturbation theory, 

>2x(2) _ ^(2) 



2\(2) 
)p ■ 



The first term is simply 



(fi 2 )i 2) = tt^ (s 2 , Df\s^ = M 22 + M k 22 (k - u h ) 



with 



and 



M, 



22 



\ S 2\ 



S2, 



,(2) 



*2 



M. 



22 



P2 



«2, 



(2) 

-2o;2 + 2^03 



s 2 



(36) 

(37) 

(38) 
(39) 



To calculate the other contributions, let us first introduce the following notations (SA 
component dropped - it is decoupled in all orders; see also Appendix iBl) 



(40) 



In Eq. (|40l) . ip 2 is the unstable mode of the embedded ANO vortex (see Eq. ([34]) ) and ipl is its 
conjugate mode (it can be seen easily, that the equation of s 2 is the same as that of s 2 with the 
transformation k — > — k). tps and ipA are the symmetric and antisymmetric bound modes of 
the ANO vortex, respectively, while a 3 is another bound mode of the ANO vortex, satisfying 
the same equation as ghosts. 

The dispersion relation of these modes (in the same order as in Eq. (1401 ) is given as 



1p2 


= (0,0,s 2 ,0,0,0, Of 
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= ( s o > s o > 0' 0' a o > a o > 
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-<tf,0) T 


a 3 


= (0,0,0,0,0,0,a 3 ) T 





^2 


= \j + (k-u) 2 




= \l + (k + uj) 2 


n 2 


= \l+k\ 


nl 


= X 2 A + k 2 , 


n 2 * 





(41) 



with Aj being constants (lowest eigenvalues). For second order perturbation theory, the following 
matrix elements are needed, 



M 22 , 
M 2S 
M 2A 
M 23 



1 



l S 2| 



s 2 ,U, 



1 



I S 2|||FS| 
1 

MIIIV^I 
1 

l s 2|||KH 



S2, (V W 







V^)s s ) + [s 2 ,(A^+A' 2 



(i) 



(l)^ 



(i) 
2 



A' 



(42) 



[s 2 ,B { 2 ] a 3 
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The part of the 0(e ) correction to the eigenvalue due to perturbation theory can be written 

as 



1 )p A| - 2uo 2 + 2kcu h \\- 2u 2 + 2ku b \ 2 a3 - 2cu 2 + 2ku h K ' 

The contribution to the perturbation functions will be evaluated similarly using the formulas 
well known from elementary quantum mechanics in Subsection 12.31 

Looking at the potentials in the equations of ips an d ipA ( see Figure [TOl in Appendix [B]), one 
can easily understand why the mode ips becomes quasi-bound (instead of bound) for values of 
(5 above 1.5 (see Appendix iBl) . The contribution of these modes to PT has to be taken into 
account when calculating the perturbation (1431) of the eigenvalue. This can be done as follows: 
the s-component of the quasi-bound mode varies slowly with the energy in a given interval, 
while the a-component varies rapidly. This makes it possible, to smoothen the contribution by 
using a low-lying quasi-bound mode, with the a-component taken to be zero and reading off 
the integral of the state density over the energy interval, in which the s-component is nearly 
constant, from the perturbation functions, and using this correction factor a to calculate the 
energy correction: 

= a]M2s]2 + ^ + ^ (44) 

1 )p A| - 2cu 2 + 2ktu h \\- 2cu 2 + 2kco b \ 2 a3 - 2cu 2 + 2kto h 1 ; 

with the matrix elementM 2 s calculated from the smoothed ips, and a is obtained from a 
comparison of the perturbative and exact eigenmodes. 

The eigenvalue fl 2 at k = u; b as a function of e (near ui = oj^) can easily be calculated 
with perturbation theory as outlined in the previous section. Second order perturbation theory 
predicts 

^U = (^) (0) + e 2 (^) {2) - (45) 

Values of (VL 2 )^ and (Vt 2 )^ obtained with the perturbation theory methods of the previous 
section and exact (numerical) values are compared in Table [5] and Figure [3 Similarly to the 
case of the vortex backgrounds close to the bifurcation, there is a good agreement between the 
exact eigenvalues and the perturbative ones, for moderate values of e, i.e. approximately for 
e < 0.2. This is an indication that the instability mode of the twisted vortex can be regarded 
as a deformation of the instability mode of the embedded ANO vortex. 

2.3 Perturbative calculation of the eigenfunction corresponding to 
the instability 

In this section, we present the perturbation functions Si/, s{ _ e , S2,e, _ e , ae, a*_ e , a^e. We will 
also compare these with their first order perturbative expressions near the bifurcation point. 

Let us first note, that for vortices close to the bifurcation (i.e. u close to Wb, £< 1), the 
lowest lying negative eigenvalue is close to k — u. Let us now choose a vortex which can still 
be treated perturbatively, but it is not very close to an embedded ANO vortex, eg. to — 0.13 
for P = 1.25, and examine its instability mode for k = uj. 

In the previous subsection, we have applied perturbation theory to obtain the eigenvalues 
of the perturbation operator. Let us now apply PT to evaluate the wave function, i.e. the 
unstable mode of the twisted vortex. The unperturbed wave function is simply the instability 
mode of the embedded ANO vortex, 

^(°) =^ 2 = (0, 0,s 2 ,0,0,0, 0) T . (46) 
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Let us now write up the first order perturbations to this mode, 

M 2S , Ar M 2A 



■0 



(i) 



-N 



A| - 2u>l + 2ku] 



4s - N- 



4; 



N- 



M> 



23 



-a 3 



(47) 



A^ - 2cu 2 + 2kto h r A3 - 2ujI + 2/cw b 

where N = \\1p2W. The first order wavefunction and the exact instability mode exhibit a good 
agreement, even for not so small values of e (see Figures [8] and [U in that case e = 0.264). 




Figure 8: The twisted vortex instability mode, scalar field components for j3 = 1.25, u 
Perturbative results for si, s* and s 2 are shown with dotted lines, s 2 = to this order. 



0.13. 



Note that the perturbations depicted in Figures [EHHl correspond to a quite large value of e, 
and even so the instability mode components are still strongly dominated by s 2 . This indicates 
that the instability of the twisted vortex proceeds similarly to that of the embedded ANO 
one, i.e. the vortex starts to expand nearly homogeneously. To support this conclusion further, 
we calculate the perturbations of the current of the vortex. The components, which receive 
perturbations, are 



(48) 



where, s\ _ e = 0(e), = 0(e), and f 2 = 0(e), i.e. this perturbation does not change the zero 
component of the current in the leading order, and 



#3 



£ M i{Qt-kz) 



{(k + 2uja 3 )f 1 s 1/ — (k + 2uj(a 3 - l))f 2 s 2/ 
+ {-k + 2ua 3 )f l8 ! t - {-k + 2co(a 3 - l))hs* 2 . + 2{fl 



(49) 



/2 2 )«3/} 

where s\ _ t = 0(e), s M = O(e), a 3/ = O(e) and f 2 = 0(e) therefore = O(e) too. 

The deviation of both the eigenvalues and the eigenfunctions of twisted strings is of 0(e 2 ) 
with respect to those of the ANO ones. This suggests that the eigenmode of the ANO string 
is an energy lowering perturbation of the twisted one. The energy, Q 2 = {jp 2 , D 2 ip 2 ), of the 
embedded ANO instability mode, s 2 , on a twisted vortex background can be expanded as 

fi 2 = (ft 2 )(°) + e 2 (fi 2 )(°). (50) 
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Figure 9: The twisted vortex instability mode, vector field components for (3 = 1.25, to = 0.13. 
Perturbative results are shown with dotted lines. 



p 


M 22 


M k 22 




^max 


1.25 


0.6883 


0.2934 


-0.0188 


0.1653 


2 


1.2442 


0.5990 


-0.1088 


0.2957 


2.5 


1.4843 


0.7270 


-0.1827 


0.3396 



Table 6: Data for the calculation of the energy of the mode ijj 2 on a twisted vortex background 

Using the expansion of the perturbation operator 

D 2 = D® + , (51) 

yields 

(n 2 )^ = M 22 + M* 2 (k - u h ) . (52) 

Table [6] shows clearly the important fact that for not too large values of e, the perturbation s 2 
still lowers the energy. The e max values were calculated with k = 



2.4 Properties of the instability modes 

Let us recapitulate first, that, as shown in Subsection 12.11 the obtained unstable eigenmodes 
are all physical. An examination of the above calculated eigenvalues and the corresponding 
eigenfunctions shows clearly that the obtained instability eigenmode is a deformation of the 
instability eigenmode of the embedded ANO vortices already mentioned in Section 11.31 (See 
also Appendix [B|), in the sense that a for a given value of the z-direction wavenumber k, 
the corresponding eigenmode \& and eigenvalue are smooth functions of the twist u, which, 
for uj = co>b reproduce the eigenmodes and eigenvalues of the embedded ANO vortices. In 
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Subsection 12.31 we presented a description of this deformation based on perturbation theory, 
with the unperturbed solutions being the eigenmodes of the embedded ANO vortex. 

The most important result is that for each k G [— k m , k m ] there is one eigenmode correspond- 
ing to a negative eigenvalue (here k m , of course, depends on u). In this section we would like to 
summarize the properties of these eigenmodes, and speculate about the physical consequences 
of them. 

Firstly, let us briefly analyze the properties of the instability modes. Let us point out that 
even for values of to quite far from the bifurcation u h (e ~ 0.2... 0.3), the unstable mode 
is dominated by the S2 component (see also Section 12.31 and Figures [HHH] therein), and the 
lowest lying eigenvalue is close to k = u. If k = u then the dominant S2 component of the 
unstable mode is ^-independent. In the previous subsection we have shown, that for u « c<j b 
(see Table [6] for critical values of e at which the energy of this perturbation becomes positive) 
the unstable mode of the ANO vortex is also an energy lowering perturbation for the twisted 
one, albeit not an eigenmode. In this way we can explicitly construct a z-independent energy- 
lowering perturbation of the twisted string, which is an indication that homogeneous expansion 
of twisted strings close to the bifurcation is possible. 

The deviation from fc m ; n = u> occurs at small values of u, far from the bifurcation, on = u^, 
i.e. the ^-dependence of the instability eigenmode becomes important when the vortices have 
already quite expanded in the lower scalar field component. In our opinion, this is another 
indication that the instability corresponds to an expansion of the string. 

In order to see how the string can expand in case a completely homogeneous expansion is 
not possible, let us take a perturbation of the form of a wave packet, centered around z = 
in coordinate space and around k = fc mm in momentum space. Then, looking at the Fourier 
transform of such a wave packet one gets AxAk = 1 for the width of the wave packet in 
coordinate and momentum space, respectively. Let the amplitude of the packet be A. Then, 
the energy of this packet can be calculated as 

^ = N 2 Ki n + W] + - 5 

which can be negative if Ak is sufficiently small. However, for the wave packet to make the 
current grow (locally make the vortex resemble a vortex of lower u, i.e. let the vortex locally 
get diluted), one needs 

Ax < A 



k 



mm 



which gives 

k 



Afc~ — 



2tt 

Substituting this into the expression for the energy, at least close to the bifurcation point, one 



can see that 



£^| 2 (^-lKm<0. 



We have also checked, using the values of f2^ in and f2| obtained above, that for j3 = 2 for all 
values of oj available for our numerical methods the energy contribution of the lump can be 
negative, which suggests that the instability of the vortex corresponds to an expansion insta- 
bility, similarly to the case of the embedded ANO vortices. This result makes it plausible that 
a local (lump-like) spreading of the vortex resulting in a configuration with diluted magnetic 
flux is possible. Forthermore, a (local) expansion of the vortex can form out of a general initial 
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perturbation of the string, because the eigenmodes which are almost ^-independent are the 
ones which are close to mode corresponding to the largest negative eigenvalue, and therefore 
they are the ones that grow most rapidly in time, and a general initial perturbation (like the 
lump discussed above) overlaps with these modes. 

The existence of the ^-independent energy-lowering perturbations for strings close to the 
bifurcation and local lumps for all strings independently of u clearly show that at least at 
the outset of the instability an expansion of the string is possible. In the case of embedded 
ANO vortices it is known [3J [TQj , that they are unstable against expansion. The similarity of 
the eigenmodes of twisted and embedded ANO strings supports that at least initially, twisted 
strings also start to expand. However, in the linear approximation, the perturbations of the 
current components, Eq. (1481) and (1491) . are harmonic in z, and thus the conserved current of 
the twisted string remains localized. This changes the dynamics of the expansion, compared 
to that of the embedded ANO strings considerably; the expansion cannot go on indefinitely 
until the current is carried away by radiation. Radiation effects only occur in higher orders 
of perturbation theory, and thus the framework of the present paper, i.e. linearization, is not 
sufficient to draw definite conclusions about the long term dynamics of the instability, and the 
resulting final state. 

The stability of the global version of current carrying strings, in a theory with two scalar 
fields, has been analyzed in Ref. [27]. It has been found that the string is stable against 
homogeneous and only ^-dependent perturbations, and above a critical value of the (global) 
current unstable against time dependent ones. For twisted vortices in the semilocal model, it 
is the magnetic flux in the string that is repsonsible for the spreading instability. 

Ref. [H] suggests another scenario for the instability of semilocal twisted strings, based 
on the harmonic z-dependence of the S , f/(2)-current flowing in the string (see Eqns. fl48l) and 
fl49l) in the previous subsection), and the fact that a larger current corresponds to a more 
expanded (lower uj) string. There, an analogy with the Plateau-Rayleigh instability of a fluid 
stream (droplet formation) is emphasized. The analogue of the Plateau-Rayleigh instability, as 
suggested by Ref. [16] is the break-up of the string into small, droplet-like string segments. 

As mentioned above, there are good reasons to believe thet for not too large values of the 
current, twisted strings start to expand initially. For large values of the current where the 
nearly homogeneous expansion is not an energy lowering perturbation, the width of the string 
is already quite large, therefore, in neither case can one conclude that the string will initially 
start to break up. 

Recently, in Ref. [28] the dynamics of twisted vortex loops (vortons) has been studied in 
a model with two scalar fields. There, it has been found that such a loop is stable against 
axially symmetric perturbations (eg. expansion), and unstable against non-axially symmetric 
ones. Longer time evolution shows that the loop breaks up into smaller pieces. In our opinion 
however, this result in itself is not decisive to support the breakup scenario of twisted strings 
in the semilocal model. Besides the different initial dynamics, anoter reason why breakup of 
the string is unlikely is that, in order for the string to break up, string ends and droplet-like 
configurations must form. Clearly if a string breaks, its magnetic flux must end in something 
like a magnetic monopole resulting in a large energy contribution. We expect the situation to be 
rather different when non-Abelian gauge fields are present, since then a possible configuration 
for a string piece would be the famous dumbbell of Nambu [T8] . 

Another interesting question examined in Ref. [E] is the stability of solutions with periodic 
boundary conditions at z = 0, L (this is an approximation of a vortex loop, if L is large). In 
the semilocal model, periodicity constrains the twist, so that uL has to be an integer multiple 
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of 2vr, 

LuL = 2np, peZ. (53) 

In this case, for not very small values of u (i.e. not very far from the bifurcation), the k = uo 
mode corresponds to a negative eigenvalue, therefore these short strings with periodic boundary 
conditions are unstable in the semilocal model. 

An interesting observation of Ref. [16j is, that considering twisted strings as solutions of the 
SU{2) x U{1) (electroweak) gauge theory, the instability found in the semilocal model can be 
excluded by imposing periodic boundary conditions. In the electroweak model, with a suitable 
gauge transformation involving the electromagnetic and the third SU{2) component vector 
fields, the ^-dependence of the string can be gauged away and therefore periodic solutions of 
arbitrary length exist. However such periodic strings are not gauge equivalent anymore to a 
piece of a twisted string, the latter being a discontinuous configuration with infinite energy. 
Looking at the gauge transformation connecting the twisted and the untwisted string, the 
generator of this transformation is proportional to z. When z is made periodic, it ceases to be 
a univalued function, and thus, the vector field component proportional to dz ceases to be a 
pure gauge. 

The gauge transformation which removes the twist of the vortex in the gauged SU(2) x 
U(l) model renders the homogeneously expanding mode ^-dependent, and thus, this perturbed 
solution is excluded by the periodicity condition. This means that these small periodic strings 
are stable with respect to purely U(l) perturbations, however, perturbations in the SU{2) 
non-Abelian sector can still cause instabilities. Whether such instabilities exist, could only 
be decided considering all perturbations in the SU{2) x U{1) gauge theory, as in the case of 
electroweak strings [21]. Ref. [21] performs the linearized stability analysis of electroweak Z- 
strings. There, Z-stings are found to be unstable unless (5 is small or 9 W is large, and the 
instability mode is ^-independent, occurs in the non-Abelian sector, and is decoupled from the 
other field components. However, these properties of the perturbation equations cannot be 
generalized to the case of the two-component strings of Ref. [16j. 

3 Conclusion 

We have found that twisted strings in the SU(2) g \ ohal x U{1) extended Abelian Higgs (semilocal) 
model possess a family of unstable modes, parametrized by their wave number in the z direction, 
k G [— k m , k m ]. This result is in agreement with that of Ref. [IE]. We have also obtained a semi- 
analytical description of twisted strings based on a suitable perturbative expansion near their 
bifurcation point with the embedded ANO solutions. The perturbative solution provides a 
surprisingly good quantitative approximation for twisted strings themselves, as well as for the 
instability eigenvalues and eigenfunctions. This way the instability modes of twisted strings 
are obtained as deformations of those of the embedded ANO solutions. 

Based on the similarity of the unstable modes of twisted and embedded ANO strings, we 
suggest that the initial time dynamics is also similar, i.e. at the outset of the instability twisted 
strings start to expand (nearly) homogeneously, just like embedded ANO ones. This idea is 
also supported by the fact, that the (homogeneous) lowest energy eigenmode of the embedded 
ANO string is also an energy lowering perturbation for twisted strings with u not far from the 
bifurcation u h , although not an eigenmode. For values of the twist uj <C u h , the formation of a 
growing lump leading to a local expansion of the string, is also an energy lowering perturbation. 
This suggests that for these strings (where the condensate in the core is already quite sizable), 
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a local expansion is thus a plausible scenario for the initial dynamics of the instability. Twisted 
strings carry a global current. Although the initially homogeneous distribution of this current 
is changed into a ^-dependent one, the current remains localized in leading order perturbation 
theory. Therefore the expansion of the string cannot proceed indefinitely, unless higher order 
radiation effects do indeed remove the current from the string. 
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A Details of the perturbation equations 

In this section we present some of the omitted details of the calculations. Although in other 
parts of the paper we have only studied the n — 1, m — case, here we present the general, 
arbitrary n and m case. 

A.l The perturbation operator of the twisted vortices 

Let us apply the Fourier transform (1291) to the equations (1261) . A purely harmonic mode can 
be written in detail as 



<tyi(z, t; Xi) = e i(m " fez) 50i(A;, ft; Xi ) 
5(j) 2 (z, *; Xi) = e t(m - (fc -" )z) 50 2 (£;, ft; Xi ) 
5A fl (z,t;x,) = e lint - kz) 5A fl (k, ft; xA 



S(f)l(z,t; Xi 



i(Qt—kz) 



)j ( k, ft; Xi 



%(z, t; x t ) = e l{nt - {k+uj)z) 5<j)* 2 {-k, -ft; a*) (54) 



with the index i running over 1,2. The variables are real functions, therefore 

A^[k, ft, Xi) = A (—k, —ft, Xi). 



(55) 



Substituting these into equations (1261) yields the perturbation operator M. of equation (1301) : 



M 



(v, 




Vi 




A x 


Bx 


0\ 




V\ 


V[* 


vi 


A\ 


B\ 





v 2 


v 2 


v 2 


u 2 


A 2 


B 2 





v 2 * 


v; 


w 2 


vi 


A* 


B* 





At 


Ai 


A* 


A 2 


v 3 










Bx 


B\ 


B 2 





-Ds 





v° 

















v 3 J 



(56) 



with 



v 2 
v; 



k 2 + V s + 2iAidi + 2kA 3 + Wx , 

k 2 + V S - 2iAidi - 2kA 3 + Wx , 

(k - oo) 2 + Vj + 2iAidi - 2(oo - k)A 3 + W 2 

(k + uof + V f - 2iAidi - 2{u + k)A 3 + W 2 

k 2 -d 2 + 2\6\ 2 , 
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where 



and 



-d\ + A* + Al + /3(|0| 2 - 1) W* = (/3 + 1)0*0* 





= C9 - 1)0! 


w 2 




Vi 


= (/3 + 1)010; 


v 2 


= (/? + l)020t 




= (/? - 1)0X02 


v 2 


= 09 - 1)020! 


Ai 


= 2Ai0! + 2z9i0l 




= 2A 4 2 + 2z^0 2 




= 2A 3 0i 


i5 2 


= 2A 3 2 



The Fourier-transform takes the equation (1281) of the ghost modes into 

while a gauge transformation takes the form 

SAi -> 5Ai + diX, 5A 3 -> 5A 3 - ifcx 
It is useful to introduce complex coordinates 



(57) 



(58) 



^-=(A r + -A* 
V2 r 



Fourier expansion in the angle variable (in cylindrical coordinates x 1 = r. < 
omitting the sum over £ yields 



(59) 

2 - 0, see Eq. (ETJ)), 



50! (k, SI) 


= s^+W , 


m 


-k, 


-SI) 


— b X,-£ e 


5(fi 2 (k, SI) 




m 


-k, 


-SI) 


— e* ^-*( m -0^ 


5A + (k, SI) 




<L4_( 


-k, 


-SI) 


= _ ia ^ e <«+i)* 


6A (k,Sl) 




* / 

a o,-A 


-K 


-SI) 


= a ,i(k, SI) , 


5A 3 (k, SI) 


— n P m 


a t-e( 


-k, 


-SI) 


= a 3/ (k, SI) . 



(60) 



Substituting this into the equations of motion, (|26j) assumes the form of an eigenvalue problem 
fl32l) with the operator 



Me 



fD 1 


lh 


V 


V 


Ax 


Ax 


Bx 


o\ 


Ut 


D\ 


V 


V 


A[ 


Ax 


Bx 





V 


V 


D 2 


u 2 


A 2 


A 2 


B 2 





V 


V 


U 2 


D* 


A 2 


A 2 


B 2 





A, 


A 1 


A 2 


A> 2 


D 3 











A[ 


A x 


A' 2 


A 2 





D* 








Bi 


B x 


B 2 


B 2 








D 4 





\o 




















D 4 ) 



(61) 
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with 



D l = D s + 
D{ = D s + 
D 2 = D S + 
D* 2 = D S + 
D 3 = D a + 
D* 3 = D a + 
D 4 = D a + 



n{\ — a) + 



n(l — a) 



m — na + 



+ k 2 + 2ktoa 3 + Wi 



+ k 2 - 2kua 3 + Wt 



T 

m — na 



+ (k- uo) 2 + 2(k - uo)uoa 3 + W 2 
+ (k + uj) 2 - 2(k + u)ua 3 + W 2 



-If 



with 



D s = -V r 2 + ^ + /5(|/| 2 -l) 

W t = (P + l)f 2 , 

D a = -V 2 r + k 2 + 2\f\ 2 



and 



Ul 


= (/?- i)A 2 


c/ 2 


V 


= (/5+l)/l/2 




A 1 


= -V2(n-fa-a)) 


A 2 




= ^(£ + ^(l-a)) 




By 




5 2 



C/ 2 = {fi - l)f 2 
V = (J3- l)/i/ 2 



m — na 



r 

m — na 



B 2 = 2oo(a 3 - l)f 2 . 

The expansion of the gauge transformation generator function can be chosen as 

ifit 

X = Xtz ■ 

Using this expansion, the ghost mode equation (1571) assumes the form 
Gauge transformations satisfying the above equation act on the fields as 

Sal —> S a 



h 



iXefa 
i 



a e - 



V2 

0-31 -> a 3,i ~ ikxe 



Xe + 



IXlfa 

i 



a_ e — > a_ e + ^_ 
a>oi -> a>z± + i^Xe ■ 



Xe 
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A. 2 The (3 — ► oo limiting case 

In the (3 — > oo case the linearized form of the constraint 0*0 a = 1 is 

W. + W; = 0. (66) 
This can be taken into account by the substitution 

si,e — tecosO + tij, s{ _ e = t e cos6 - t^ e , 

(67) 

S2,e = —U sin 6 + t 2) e, s 2 _ t = -t e sin 6 - £ 2 ,i, 

and the equations of tf, and t 2j ^ can be obtained with a somewhat lengthy but straightforward 
calculation. Here we present only the t = case: 

-«)' = + Ao*o + £*o*i + Fto*2 - sin(2fl)(a + oj) + sin(2#)a 3 ,o (68) 

with A t0 = —(m—n)(m+n—2na) cos(2#), Ao = k 2 — Q 2 — u 2 (l — 2a 3 ) cos(20), Ao = 2£xua 3 cos6> 
and F t0 = 2kui(l — a 3 ) sin#. Similarly 

-K)' = + CtA + Ai^o + Ai*2 - V2 cos 66' (oq - a*) (69) 

with = -| [m 2 - ri 2 - 4(m - n)nacos 2 6 + (m 2 - n 2 ) cos(20)], C n = 1 + k 2 - f2 2 - ^ - 
cos(2#) + 2a; 2 a 3 cos 2 # - ^cos(2fl) - (6'f , Ai = 2kua 3 cos6 and E tl = sin(20), and 

-K)' = ~r*2 + C a t 2 + D a t + 3b*i + v^sin^'(ao - a*) (70) 

with A t2 = -\ [n 2 - m 2 + 4(m - n)na cos 2 6 + (m 2 - n 2 ) cos(20)], C t2 = 1 + k 2 - Vl 2 + + 
cos(20) — 2co> 2 a 3 sin - ^ cos(20) - (0') 2 , A 2 = -2ku(a 3 - 1) sin0 and E t2 = sin(20). 

A. 3 The expansion of the perturbation operator in the bifurcation 
parameter 

Let us first expand the components of the perturbation matrix in Eq. (|6~TT) to second order in e. 
Here we present only the terms used later on in this paper. The differential operators expand 
as 

= -VI + ( "r +f)2 + (* - u h ) 2 + (3(Uf y ) 2 - 1) , 
D?=0, 

= 2u h u 2 - 2ku 2 - 2(m ~^ +£)na(2) + 2u h (k - . b )4 2) + /f + (2/3 + 1)(/ 2 (1) ) 2 , 

(71) 



D r = -v 2 + (m -^ 0) ' £)2 + (* + - b ) 2 + /^((/i (0) ) 2 - D , 

A* (1) = o, 



(2) _ o,.„ ^ oiv... _ 2(m-na(°)-l)nq( 2 ) _ Q _ ^ ^ _ ^ (2) ^ Wo) /(2) ^ ^ ^ 1 v ^ (1))2 

(72) 



D?*> = 2iu h LU 2 + 2ku 2 - ^ '- 2cu h (k + cu h )a^ + 2/3f^f^ + (2(3 + l)(/ 2 l 
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Scalar-scalar interaction terms: 

V® = 0, 
v m = , 

Scalar-vector interaction terms: 

while Af = A' 2 W = Af = A' 2 (2) 

This e-expansion of the perturbation matrix elements can be used to obtain the unstable eigen- 
value of twisted vortices close to the bifurcation point via perturbation series. 

B Perturbations of the ANO vortices 

The equations describing local vortices can be obtained from the twisted vortex equations ((9)) 
by taking / 2 = a 3 = 0. With the same numerical methods as used in the case of the twisted 
vortices, local vortex background data can be obtained, see Table Here, the and a^ 2 ' 
is defined by the behavior of the vortex profile functions at the origin, / ~ /«r + 0(r 2 ) and 
a ~ oP^r 2 + 0(r 4 ) (shooting parameters). 









1.25 

2 

2.5 


0.92418 
1.09935 
1.19677 


0.53485 
0.61657 
0.65959 



Table 7: Local vortex background data 





Symm. mode 


Anti. mode 


a 3 mode 


p 




X 2 


X 2 

A a3 


1.25 


1.82404 


1.62442 


1.62442 


2 




1.76100 


1.76100 


2.5 




1.81813 


1.81813 



Table 8: Parameters for the £ = ANO vortex perturbation modes for (3 = 1.25, 2 and 2.5 

The perturbations of local vortices can be described using equation flBTl) with uj = and 
f2 = o>3 = 0. The ctse component decouples, and it's equation becomes that of the ghost modes, 
Eq. fl64l) . The s^i and s% _ £ components also decouple; their equation (fT4l) has already been 



^ = (/5 + i)/ 1 (0) / 2 (1) , W = o, 



v m = {p- i)/r/ 2 i 



(o) ,(i) 



V 



(2) 



0. 



(73) 



m — na^ 
r 



f (i) 

32 



B. 



(o) 



(2) 



0, and 



(74) 



5. 



(i) 



(i) 



(75) 



26 




S mode scalar 
AS mode scalar 
Vector 




1 ' 1 1 1 ' 1 

2 4 6 8 10 12 

Figure 10: Potentials in the ANO vortex perturbation equations for (3 = 1.25 (left), (3 = 2 
(right) and (3 = 2.5 (below) 

given in Section [TJ The instabilities arising from the extension of the model with a second 
Higgs field are described by these equations. For further details see [31 [9]. In the £ = case, 
there is a further decoupling in the (si,0) s* , a , a^) sector: symmetric: ip$ — ( s i) s i) a o ; a o) an d 
antisymmetric if) a = (si,—si,ao,—ao) modes are decoupled. For values of (3 < 1.5 there is a 
bound mode in each of these two sectors. For (3 > 1.5 the symmetric mode becomes quasi- 
bound. The parameters for these modes (normed in a way that s[(0) = 1) are summarized in 
Tabled There, A^ Sa3 denotes the lowest eigenvalue (corresponding to k = 0). Figure [TOl shows 
the potential terms for So and ao in the symmetric and antisymmetric sectors, respectively, for 
(3 = 1.25, 2 and 2.5. These explain it clearly why the mode ifts becomes quasi-bound for higher 
(3 (i.e. for (3 > 1.5, there is no bound mode in the symmetric sector, but for 2 < Q 2 < 2(3, the 
scalar channel of the scattering problem is closed). 

The gauge dependence of the eigenmodes and the relation of ghost modes to these is rather 
straightforward in the case of the ANO vortices. In the analysis of the spectrum of these 
vortices, in Ref. [H] it has been noted that if one considers the perturbation problem of the 
vortices in 2+1D (i.e. ^-independent perturbations), then the ghosts cancel the perturbations 
of A , and some modes of the gauge fixed perturbation operators. If we allow ^-dependent 
perturbations, instead of SA a combination of 5 A and 5A 3 is canceled. Ref. [TT] also identifies 
the non-physical mode with the antisymmetric mode in the £ = case, based on the small 
difference in their numerical eigenvalues. Here, we would like to supplement this identification 
of the canceled mode with some analytical calculation. Examining Eqns. fl65"]) . one can easily 
conclude that in the symmetric case, the gauge generator drops out of s± = (s^o + s* )/2 and 
^0 = ( a o + a o)/2- On the contrary, such a cancellation does not happen in the antisymmetric 
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case, thus the ghost eigenvalues shall appear in the spectrum of the antisymmetric mode. 
On a twisted vortex background, the couplings are more complicated, and therefore such a 
straightforward analysis is not possible. 
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